Coupling DuMuX and DUNE-PDELab to investigate evaporation at the interface between Darcy and Navier-Stokes flow

An implementation of a coupled Navier-Stokes/Darcy model based on different Dune discretization modules is presented. The Darcy model is taken from DuMux, the Navier-Stokes model is implemented on top of Dune-PDELab, and the coupling is done with help of Dune-MultiDomain together with some project-specific auxiliary code. The Navier-Stokes model features one fluid phase, the Darcy model two fluid phases. Each fluid phase may be composed of two components, in addition, non-isothermal processes are considered. The coupling between free and porous-medium flow uses a sharp interface between both subdomains and conserves mass, momentum, and energy by accounting for the corresponding fluxes across the interface. A cell-centered finite volume method (FVM) is combined with a marker and cell (MAC) scheme. It solves the coupled problem in one monolithic system using a Newton method and a direct linear solver. Numerical results demonstrate the basic functioning and a lab-scale reference application.